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ABSTRACT 

This paper deals with predator-prey model having Holling type III functional response. The prey population is 
stage structured consisting of immature and mature stages and the predator population is influenced by the resource 
biomass. Dynamical behaviors such as positivity, boundedness, stability, bifurcation of the model are studied analytically 
using theory of differential equations. Computer simulations are carried out to prove the analytical result. It is noted that 
influence of resource biomass on the predator population may lead to the extinction of predator at a lesser value of maturity 
time in comparison to the absence of resource biomass. 
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1. INTRODUCTION 

One of the most powerful tools for understanding the predator prey relationship in ecology is to understand the 
dynamic relationship between predator and their prey. Many researchers have studied the system in depth [1, 7, 14, 17, 18, 
19, and 22]. The Stage structure population where the individual member have a life history that takes them through two 
stages, immature and mature have received much attention in mathematical modeling of ecological system. In all these 
studies, the maturity age is represented by time delay thus the resulting in the system of retarded functional differential 
equations. Zang ET. ah, [26] in their paper investigated the behavior of predator -prey population with prey as structured 
population. Song and Chen [20] studied a two species competitive stage structured population model with harvesting for 
prey. They obtained stability conditions and threshold of harvest effort for population survival. 

It may be pointed out here, that most of the above studies are based on the traditional predator-prey models with 
either prey or predator stage-structured or both. In nature, there are many cases where predator population dynamics is 
influenced by the presence of an additional resource (which may or may not be a secondary prey). This can be illustrated 
with the example: stone martins, relative of the weasel, are extremely fierce and dangerous predators, and often take prey 
like squirrels, pike, voles, hares, etc. [2]. They live in forest, especially evergreen ones and spend much of their time up on 
trees, jumping from one place to another, climbing up and down and rarely reaching the ground. They build a den in an 
abandoned hole in a tree (resource biomass) and hence produce an adverse effect on the growth of the trees. This example 
shows a relationship between resource biomass and predator population, although the predator feeds on prey only. 
Consequences of moose acting as additional prey in a Wolf-Caribou system have been described by Bergerud [4] and 
Bergerud et al. [5]. Freedman and Shukla [9] analyzed a predator-prey system where the resource dynamics affects the 
predator-prey system. Freedman et al. [8] analyzed a ratio-dependent predator-prey model where the predator population is 
influenced by the presence of a resource. 
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Waryano Sunaryo, et. al., [21] studied an ecological model with a tri-trophic food chain composed of a classical 
Lotka-Volterra functional response for prey and predator, and a Holling type-III functional response for predator and 
super-predator. Agarwal and Pathak [3] studied the effect of harvesting on dynamics of prey predator model with holling 
type III functional response. 

It may be noted from the above investigations that behavior of resource biomass on predator- prey stage structure 
model with Holling type III functional response is not studied yet. Therefore in this paper our focus is to model the effect 
of resource biomass on predator prey population with Holling type III functional response where prey population is stage 
structured. 

2. MATHEMATICAL MODEL 

In this paper Holling type III predator functional response on predator - prey - resource model with stage - 
structure for prey as the system of following differential equations: 

Xi (t) = ax m (; t ) - yx t (; t ) - ae~ n x m ( t-r ) , 


x m ( t ) = ae~ yr x m ( t-T ) - (3x1 ( f ) “ 


cx;„(t)y(t) 

x i(Oe l +1 


( 2 . 1 ) 


m = y(t ) 


-d(r{t)) + 


Kit) \ 

l + e i x i(t) J 


m = r(t) 


1- 


r(t) 


L(y(t)) 


x m (t) = </> m (t) > 0, - T < t < 0 and x ; (0) > 0, y(0) > 0, r( 0) > 0, 


Where x r (t) and x (t) are the densities of immature and mature prey populations respectively, y(i ) is the density 
of the predators and r(t) is the density of the resource biomass. 

In mathematical model (2.1) all the parameters are positive under the following assumptions. 

HI. OC is the proportionality constant for immature prey population to mature prey population, T is delay period 
at which immature population transfer to mature prey population, y is death rate of immature prey population and (3 be 
the intra specific interaction. 

The term CXe n x iH (t — T) describes that immature prey population born at time (t— f ) and surviving at the time 

t and further covert immature to mature prey population, c , e, and / are constants and are defined as the capturing rate 
of predator, efficiency rate of predator on prey and conversion rate of predator, respectively. 

H2. d (r) = d 0 - d , + d { e~ r Is the death rate of predator, which decreases as the density of resource biomass 
increases, 

d( 0) = d 0 > 0, d\r ) <0, r > 0 ■ 
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L(y) = L 0 - L, + L t e~ y , Is the carrying capacity of the resource and it decreases as y increases. L(0) is the 
carrying capacity of the resource in the absence of predator y. 

Where, 

L{ 0) = L 0 > 0, L, e (0, L 0 ), And 
lim L(y) = (L 0 - L t ) = (say) > 0. 

y->°o 

Now for continuity of initial conditions, we require 

0 

x l (0)= [ ae ys </> m (s)ds, 

( 2 . 2 ) 

The total surviving immature population from the observed births on - r < t < 0 . 

With the help of equation (2.2), the solution of the first equation of system (2.1) can be written in terms of 
solution for x (t) as: 


t 

x j (t)= j ae~ /lr ~ s) x m (s)ds . 

M (2.3) 

Equations (2.2) and (2.3) suggest that, mathematically no information on the past history of x ( t) is needed for 
the system (2.1), because the properties of x i ( t ) can be obtained from (2.2) and (2.3) if we know the properties of x ( t ). 


Therefore we need only to consider the following system of equations, 

cxl (t)y(t) 


x m (t) = ae yr x m ( t-T )- /5x m (t) - 


xl(t) e] +1 


y(t)= y(t) 




fx 2 m (t) y 

1 + ( t ) J 


(2.4) 


r(t) = r(t) 


1 - 


r(t) "j > 
L(y) J 


x m (t) = <t> m (t) > 0, - T < t < 0, And 3,(0) > 0, r(0) > 0. 

3. BOUNDEDNESS OF SOLUTIONS 

Lemma 3.1 


f (Y,P ^ XV /y ^ p 

The set — J ( x y r) : 0 < x < — — , 0 < — — + — < , 0 < r < L > attract all solutions initiating in 

l m P c f c/58 ’ °J 

the interior of the positive region. Where 8 = min [ae n , d 0 - t/, } ■ 

Proof: 


Consider the following equation. 
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i m (r) = ae~ rt x„(t-T)~ fix m (t) 


We obtain that. 


2 /, \ cx m (t)y{t) 


+1 


X m(0^ 


ae 


-yr 


P 


Now form fist and second equations of system (2.4), we obtain the following 


clt 


ik+Z 

c f 




ae ^(t-T) j. 3x m 2 y(d 0 -d x ) 


f 


a 2 e~ 2n J3x w 2 y(d 0 -d l ) 


pc 


a 2 e ~ 2 ^ 


-- 


/ 


yido-d,) 


Pc c f 

- min {«<?~ rr , (<i 0 - d t )} 


aV 2 " 


pc 


f x y ^ 

— + — 


/ 


2 —2 YT 

This implies that 0 < Z_ + Z< Z_£ where S = min {ae ~ n , (J 0 - </,)} . 

c / 

Now from the third equation of system (2.4), we obtain 


r(t) < r(t) 


V L o J 


This implies that, lim sup r(t) < L 0 . 

f-> 00 

This completes the proof of the lemma. 

4. EQUILIBRIUM POINTS AND STABILITY ANALYSIS: 

For this system there exist only five positive equilibrium points which are given as: 
E 0 (0, 0,0) , E t (x ml ,0,0) , E 2 { x m2 , y 2 , 0) , E 3 ( x m3 , 0, r 3 ) and E*(x m ,y\ r ) . 

Where, 


ae 


ae 


and r 3 = L 0 . 


p '" 3 p 

The existence of equilibrium points £ 0 (0,0,0), E x (x ml , 0, 0 ) and E 3 (x m3 ,0,r. 
The existence of point E n (x , , y , , 0) is given by the equations. 


) are obvious. 
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ae~ yt - Px ml (t)- 


x l, 2 (t)e l + 1 


0 , — d q + 


fx 2 m2 (t) 
l + e 1 x 2 m2 (f) 


0. 


From these equations we obtain that the point e 2 (x 2 , y , 0) exist if following conditions hold, f > a 


ae 


- yt 


> c(rj 2 + 1 )tj , where ^ = 


P V ' ’ ' ' f ~ d 0 e 1 

Now equilibrium point £*(x^,y*,r*) exist if the system of equations, 


ae n P x m (t) ' 


-d ( r(t )) + 


cx m (t)y(t) 
x 2 m (f)e, + 1 

(0 


1 + e i X m 0) 


0, 


= o, 


1-^=0, 

L(y ) 

Has a positive solution. From second and third equations of above equations we get, 
d (L( y)) 


2 

= 


= 7i{y), (say) ■ 


{f~e id (L(y))} 

Using this function in first equation, now we take, 
ck ( y ) ^ 


p(y) = 


P + 


X (y)e l + 1 


n (y )- a 2 e 2yT ’ 


then. 


(i). F (0)= fi 2 7r(0)-a 2 e- 2rr < 0, as ^(0) < 
cjc(K) A 


a 2 e- 2yi 


P 2 


and 


(ii). F(K) = 


P + 


7r(K)e l + l 


7T (K)— a 2 e lyz > 0, where K is the maximum value of y . 


Thus there will exist y such that F ( y* ) = 0 for y* e (0 , K ) ■ 

Now the sufficient condition for uniqueness of E* (x * , y * , r*) is given as, 

dF (y) 

— r Z r> 0- 

d(y) 


The Dynamic behavior of the equilibrium points can be checked by the Jacobian matrix. 


e and 

0 c 1 
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A 

0 

-yd' (r 

L( y) 

Equilibrium point E Q (0,0,0) is trivial and its characteristic roots are given by equation = CXe ,r+A>r t 
k 2 = —d 0 and k^ =1, showing equilibrium point E 0 (0, 0. 0) is unstable in direction of x m — r and stable in the 
direction of V , and hence saddle point. 

The characteristic equation for £' | (x ml , 0,0) is given as, 


ae 


~(y+A)T 


/(E) = 


■ 2 P x „ 


2cx „,y 

(1 + e 1 *, 2 „) 2 


■k 


cx„ 


2 f X n,y 
(1 + e^,;) 2 

0 


d + e, x 2 m ) 


-d (r) + 


K 


d + e.x 2 ) 

r 2 L'(y) 

t L (y )] 2 


■k 


(ae~ (r+A,T - 2p Xl -k) 


-d o + 


fx? 


(1 + e,Xj ) 


— k 


(1-1) = 0 


Then the characteristic roots of this equation is given by equation 1 = ae ( - r+A)T -2 jSx t , 

fx 2 

and 1 as— J h — >0 it is unstable in direction y — r and as Re A < 0 the eigenvalue ae 

(l + e lXl 2 ) 

so it is stable in direction X m , so E x (x ml , 0, 0) is saddle point. 


~d o + 


~(r+X)r 


f x i 2 

(l + e lXl 2 ) 
- 2j3x 1 < 0 


The characteristic equation for the point E^(x m2 , y 2 ,0) is given as, 

(l 2 + X t k + X 2 +e~ lT {X 3 k + X 4 })(k-1) = 0, 

Then the eigenvalue in the direction r is 1 So it is unstable in this direction and others eigenvalues are given by 
the following equation: 

(l 2 + X,A + X 2 +e~ lT {X 3 A + X 4 })= 0, 


Where, 


Xj — 2px m2 + 


2c x 


(1 + e lX ,„ 2 ) 


+d - 

2\2 T “0 


f X ,r 


(1 + e ! X m2 ) 


x 2 = 


2 f cx Jy 2 

(l + e lXm2 2 ) 3 


X, = -ae~ 


And 


X, = -ae' 


d o ~ 


fan 


(1 + e ! X m2 ) 


\ 
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The stability analysis of this point seen by the theorem (4.1), so this point is saddle point which is stable only in 
the direction X m2 ~y 2 - 


Now for the point E , (x m3 ,0, r, ) characteristic equation is given as. 


(ae~ lr * X)T - 2px x -A) 


-d 0 + 


fx i 2 


(1 + ejXj) 


A\(l-A) = 0 


fx 2 

Then the eigenvalue in the direction of y is given as — d H - — — > 0 is positive and in the direction of 

(l + e lXm3 2 ) 

X m3~ r 3 the eigenvalues are negative, the eigenvalue corresponding to X m3 direction is given by equation 
A = ae~ ir+A)I - 2 J3x m3 , putting the value of Xm , = — it become A = ae~ n (e~ AT - 2) . 

P 

Suppose that Re A >0 then we calculate the real part of this, we get Rc/l = (Xe 71 (e7 lV ' llA cos(r 111! A) — 2) < 0 a 
contradiction. Hence Re /l < 0 and in the direction of r, eigenvalue is 1, so this point is saddle point. 


Characteristics equation for the point h (X m ,y ,T ) is given as: 

<j> (2, r ) = A 3 + B t A 2 + B 2 A + B 3 + e AT {B 4 A 2 + B S A + B 6 } = 0, 
Where B t , B 2 . B 3 , fi 4 , B s and B b given as, 

b, . Ln - 1 + JrL + 2 fi x : + 2 <i ' _ SE 

1 Uf) H - (l + e.O l + 


(4.1) 


B 2 = 2 /5x m + 


2 cx* m y‘ 


d + e.O 2 


~ d ( r *) + \ 


A 


« 2 V 


+ e x x m 


1- 


2 r 


. 


L (y') 


2/AI + 


2 cx* y * 
d + e,x: 2 ) : 


— d ( r ’)■ 


x, 


(1 + e ! X m>, 


2/c.<V , W*V(r*)L'(y*) 


(l + e^* 2 ) 3 


[A/)] 2 


B , 


d'(r*)L'(y*)y*r 


2 ( 


[i(y*)l 2 

2r* 

L(y‘) 


2 + 


2 Xy* 


(l + ei x* 2 ) 2 


2 Px" m + 


2 AA* 


(1 + e, x* 2 ) 2 


K 2 ^ 


-^(T , )+ , fX - , 2 
l c.x 

1 m J 


2r* 

L(y'). 


2 Ax* y* 


(1+ e, x* 2 ) 3 


B^-ae 71 = 


-p<- 


* * \ 

cx m y 

(1 + ejX ,; 2 ) 2 


5 5 = ae n 


-d (r* ) + 1 ~ 2r » + 

v 7 L(y ) i + c 


: 2 A 


0*1 + 




a i * 2 \ 2 

+ e i x m ) ) 


v ' i(y ) l + e, 


2 A 


and 
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B 6 =-ae n 


d'{r')L\f)yr 2 

[ Z /(/)] 2 


-P*m ~ 


cx rn y 


(1 + e^: ) J 


d'(r’)L\y*)y*r* 


[L\y*n 


To show the positive equilibrium E*(x* m , y’,r*) is locally asymptotical stable for all T > 0 , we use the following 
theorem [16]. 

Theorem 4.1: A necessary and sufficient condition for E*(x* n ,y*,r*) is locally asymptotically stable for T > 0 
is. 


(I) The real parts of all roots of 0( A, 0) = 0 are negative. 


(II) For all real b and 2">0, 0(ib, T) ^ 0 , where i = yf— 1 . 

Theorem 4.2: Then the positive equilibrium point Z?*(x^,y*,r*) for the system (2.4) is locally asymptotically 
stable providing following conditions. 


(i) e l d(r)> 


f 


e i x ,„+ 2 


(ii) e l d(r)> 1. 

Proof: Prove of this theorem is given in two steps as follow, 
Step-I. Substituting T — 0 in equation (4.1), we obtain, 

A 3 + (Bj + B 4 )A~ + (B 2 + B 5 )A + (B 3 + B 6 ) = 0 = 0 ( A , 0) , 


(4.2) (Bj + B 4 ) > 0 , (B 2 + B 5 ) > 0 and 


[ (BjB, + B 4 B 2 + B,B 5 + B 4 B 5 ) — (B, + B 6 )} > 0 ■ 


Then by Routh-Hurwitz criterion (4.2) has all roots are negative then in the absent of delay, E* (x* m , y*,r*) is 
asymptotically stable. 

Step-II. In this step we show that (f>{A,r) ^ 0 for A — ib , for real b , 

Let b = 0 then 

(ZHO, r) = B } + B 6 ^ 0 , it is hold II condition. 

If b ^ 0 then equation (4.1) can be written as: 

— ib 3 — iB { b 2 + iB^b + B, + e~ ,hT {-B 4 £> 2 + iB s b + B 6 } = 0(ib,T ) , 

Now equating real and imaginary parts of this equation we obtain the following equations, 

— ( B 4 b 2 - B 6 ) sin br + bB 5 cos br = b 2 - bB 2 (4.3) 

( B 4 b 2 - B 6 ) cos br + bB s sinir = b 2 B, - B, (4.4) 
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After squaring and adding the above equations we get the equation, 
b 6 + b 4 (-2 B 2 + B 2 -B 2 ) + b 2 (B 2 2 - 2 B t B 3 + 2 B 4 B 6 - B 2 ) + (B 2 - B 2 ) > 0 , as 


B 2 - B 2 —2B 2 =3\jBx m+ / CX " J ' +d(r)-- fX>n 


{ l+e i X i) ( 1 + e i X m) 


1+ A_ +M(r) _^S 

L( y) d + e.x 2 ) 2 


1 2 

J 2r 

[ + ^ 

4 1 


1 

5 


Ky 


If 

( 2r ,1 

/ 

+ 2 



Jl 

U(y) J 

V 


,2 A 


■+px m + cx m y 


— d(r) + 


K 


l + e^- 


B 2 - B 2 - 2£, > 0 , if e x d{ r) > ■ 


/ 


, condition (i) 


e i x m + 2 

B i - B s > 0 ,B 4 B 6 - B X B 2 > 0 then Bp - B 2 + 2 (£ 4 fi 6 - S,S 3 ) > 0 and 


yr 2 <i (r)L(r) 

[i(y )] 2 


- B 2 


2 r 

Uj) 


-2 /B d (r) x 2cX '"^ d } V \ + 2I5X 


3 A 


(1 + eiX 2 ) 2 l + e,x 2 J 


8 yr 2 d (r)L(r) 


l>(y)] 


p x m + 


c x „,y 


then B 3 2 - B 2 > 0 if f + p d (r): 

t(l+e lX 2 ) 2 


(f + £ \ X m ) 

P*l 


2 r 

~b(y) 


cx,„yd{ r) 
O + e.x 2 ) 2 


+ Pd (r)x n 


P x 


3 A 


1+ e i X m J 


1 + e 1 X - 


> 0 


after solving this equation we get £ 1 6?(r)>l, 


condition (ii) 

So B 2 - B 2 - 2B 2 > 0 , B 2 -2B.B, + 2B 4 B 6 - B 5 2 > 0 and b 3 2 - B 6 2 > 0 ■ 

Hence (j){ib, T) 0 for real b. Therefore the unique positive equilibrium E*(x* , y*,r*) is locally 
asymptotically stable for all T > 0 and the delay is harmless in this case. 

5. BIFURCATION ANALYSIS 

The characteristics equation of the system (2.4) for the equilibrium point E* (x * n , y * , r*) is given as, 


<p(A,r) = A 3 + B { A 2 + B 2 A + B 3 + (B 4 A 2 + B S A + B 6 )e" ir = 0. 

Let the eigenvalue of equation is in the form of A = a{x) + ib(T) and function off . Differentiating (4.1) with 
respect to T , we obtain the following equation 


i 

i 

^3 

f 3A 2 + 2AB, + B 2 ) 


f 3 AB 4 + B s ^ 

VdT ) 

^ A(A 3 + B l A 2 + B 2 A + B } ) J 


yA(B 4 A 2 + B,A + B b ) J 


If A, has only purely imaginary part a(T) = 0 and b(T) ^ 0 then A = ib(T) equation (5.1) become. 


bd AY 1 _ 

f Bp - 4b 2 B 2 + 3b 4 - 2Bpb 2 + 2 £,£, ^ 


( Bpb + 3B 6 B 4 - 3Bpb 2 \ 

{dr) 

{ (b 3 -bB 2 ) 2 -(b 2 B i -B 3 ) 2 J 


{ B 2 b 2 -( B 6 — B 4 b 2 ) 2 J 
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At T = X . . 


Re 


J/t 

(if 


q , it is verified by numerically then it is hold transversality conditions. 


The value of delay T = T c is given as. 


7 — ^arccosf ^ ^ ^ + B ^ (^1^6 + ^4^3 + ^2^5) +^6^3 ! j 2 k ?r _ £ = 0 , 1 , 2. 
c & 1 Bpb 2 +(B 4 b 2 + B 6 ) 2 

obtains by the equations (4.3) and (4.4). 

6. PERSISTENCE OF THE SYSTEM 


this value 


Biologically, persistence means the survival of all populations in future time. Mathematically, persistence of a 
system means that strictly positive solutions do not have omega limit points on the boundary of a non-negative cone. The 
persistence of the system (2.4) is give by following theorem. 

Theorem 6.1: Assuming thatCfe ' ' > C , the permanence of the solution of the system (2.4) is given by the 
following conditions (6.1), (6.2), (6.3). 


Proof: From first equation of the (2.4), 
dx 


dt 

dx n 

dt 

dx m 

dt 


■ = ae vr x,„ - Bx 2 m - cx m 

in / m m 


> ae~ VT x m - px 2 m - cx m 

> (ae~ VT - Px m -c)x m -. 


' x y A 

m y 


V 1+<? 1 X m J 


Then |j m i n f x > — > 0 ^ OC' e ^ c 

P 


( 6 . 1 ) 


Second equation of the system (2.4), 


■tL=- d (r)y + Jdl^ 

dt 1 + e i x m 


ae~ - c 

Let ltm sup y < 7 0 anc ^ lim inf x > = n. 

t->°° r->°° p 


Then above equation can be written as, 

^y_ > fd dl\ _ v , then we get 
. — . -> “o y ° 

dt l + ^/T 


lim inf y > 


fV 0 V i 


/ >lo ( a e ~ VT - C )P 


(1+ erf 2 )d 0 [p 1 + ei {ae~ vz - c) 2 ]d 0 


( 6 . 2 ) 
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Now third equation of system (2.4) 


dr 

dt 


r 

r 1- 

v 


L(y)J 



l- 



, after solving, we obtain liminf r > L 0 - L. 

t->oo 


7. NUMERICAL SIMULATION AND DISCUSSION 
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(6.3) 


Analytical studies always remain incomplete without numerical verification of the results. To facilitate the 
interpretation of our mathematical findings by numerical simulations, we assume 

d (r) = d 0 - d t + d t e~ r , d t e (0, d 0 ). And consider the set of parameter values as 


ar = 0.8,/? = 0.3,j'=0.1,r = 10,c = 4.5,/=10,d 0 = 0.4 ,d t = 0.3 ,L 0 = 10 ,L, =2,e l = 0.5 ■ For 
the above set of parameter values, the equilibrium £*(0.1003,0.5297,9.1788) is obtained. By using these parameters 
the following figure is given as, 




Time Time 


10 



o 


5 i i i i i i i i i 

0 50 100 150 200 250 300 350 400 450 500 

Time 


Figure 1: Variation of Prey Population, Predator Population and Resource Biomass with Respect to Time 



Figure 2: Variation of Prey Population With Respect to Resource Biomass at Different Initials Points 
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Figure 3: Variation of Prey Population with Respect to Predator Population at Different Initials Points 



Time 


Figure 4: Variation of Prey Population with Respect to Time at Different Values of Conversion Rate 



Figure 5: Variation of Predator Population with Respect to Time at Different Values of Conversion Rate 
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Figure 6: Variation of Prey Population with Respect to Time at Different Values of Capturing Rate of Predator 



Figure 7: Variation of Predator Population with Respect to Time at Different Values of Capturing Rate of Predator 



Figure 8: Variation of Prey Population with Different Values of Delay 
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Figure 9: Variation of Predator Population with Different Values of Delay 



Figure: 10. Variation of Resource Biomass with Different Values of Delay 

Figure 1 shows the behavior of the prey, predator population and resource biomass with respect to time, we obtain 
that after some time pre, predator and resource biomass being constant that means they got their equilibrium points. In 
Figure 2 discuss the behavior of prey and resource biomass with different initial points of them, here we see that it 
converge at their equilibriums point .In same manner in Figure 3 we see the variation between the prey and predator 
population at different initial points, we get that they converge toward their equilibrium point. 

Figure 4 and 5 describe the effect on the prey and predator population with different values of conversion rate of 
prey. As conversion rate increase the prey population decrease and predator population increase. 

In Figure 6 and 7 we see the variation of prey and predator population with change of capturing rate of predator, 
there is no change of equilibrium point of prey population but at initially it increase as capturing rate decrease, on the other 
hand predator population decrease as the capturing rate increase . 

Figure 8 described the behavior of the prey population with different value of delays, as value of delay increased 
the population of prey does not change at initially it decreased but in Figure 9 predator population decrease as value of 
delay increases but in figure 10 the value of resource biomass does not change with change of delays value. 

In the absent of resource biomass behavior of the system can be seen as follow by figure and table. 
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Figure 11: Variation of Prey Population with Different Values of Delay in Absence of Resource Biomass 



Figure 12: Variation of Prey Population with Different Values of Delay in Absence of Resource Biomass 

Table 1: Absence of Resource Biomass 


T 

x,„ 

y 

10 

0.2022 

0.263 

15 

0.2022 

0.1325 

20 

0.2022 

0.05354 

25 

0.2022 

0.005567 

31.79 

0.111 

0.0000 


Now from Figure 11-12 we obtained that as the value of delay increases value of prey population at initially does 
not change but as r > 25 value of prey slightly decreases, but in predator population decreases as delay value increases. 

CONCLUSIONS 

In this paper we studied a prey-predator-resource biomass mathematical model where prey species are delayed. 
Delay of prey is taken as the time of maturation period of prey, predator does not interact with immature prey. After 
numerical simulation we obtain following results, 

In this paper by numerical simulation for parameters value system is stable in certain conditions, here we see that 
as conversion rate of predator is increases prey population decreases but predator population increases in the same manner 
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see the effect of capturing rate of Predator on prey and predator population. 

Here we see the behavior of the system at values of delay, as in paper of Agarwal and Devi [2], they describe that 
as value of delay changes effect of delay on resource biomass is negligible but in this paper’s model we use the holling 
type III functional response as the present of this term we see as the value of delay change value of resource biomass also 
change. 
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